Wetting layer thickness and early evolution of epitaxially strained thin films 
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We propose a physical model which explains the existence of finite thickness wetting layers in 
epitaxially strained films. The finite wetting layer is shown to be stable due to the variation of the 
nonlinear elastic free energy with film thickness. We show that anisotropic surface tension gives 
rise to a metastable enlarged wetting layer. The perturbation amplitude needed to destabilize this 
wetting layer decreases with increasing lattice mismatch. We observe the development of faceted 
islands in unstable films. 
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Epitaxial deposition of a thin film onto a substrate in 
cases involving lattice mismatch is central in the fabri- 
cation of semiconductor and optoelectronic devices. The 
lattice mismatch between the substrate and the film gen- 
erates strain in the deposited film, which can cause film 
instability unfavorable to uniform fiat film growth. The 
strained film can relax either by the introduction of dislo- 
cations or by the formation of dislocation-free islands on 
the film surface via surface diffusion. Early film growth 
tends to occur via the second mechanism and we shall 
only consider dislocation free films. It has been observed 
experimentally ||l|,|[ that dislocation free flat films of less 
than a certain thickness (the critical wetting layer) are 
stable to surface perturbations, while thicker films are un- 
stable. The thickness of the wetting layer is substance de- 
pendent and decreases with increasing lattice mismatch 
strain g], e = (a^ — af)/af, where Ug and a/ are the 
substrate and film lattice constants. Above the critical 
wetting layer 3D dislocation-free islands form. Predic- 
tion and control of wetting layer thickness and an under- 
standing of early thin film evolution are important for 
the improved fabrication of semiconductor devices. 

Despite considerable efforts (see, e.g., [^ ||), the 
physics of the critical wetting layer is poorly understood, 
and the purpose of this Letter is to study its properties by 
considering the following two important questions: First, 
why is there a critical, stable wetting layer and what con- 
trols its thickness? Second, since in most cases heteroepi- 
taxial growth is done below the roughening transition, 
how does anisotropic surface tension affect the thickness 
of the critical wetting layer? Here we attempt to answer 
these questions and to study early film evolution with- 
out deposition. A later paper will look at the effects of 
deposition and long-term growth. 

We studied an elastically isotropic system under plane 
strain ||], which causes the system to be effectively two 
dimensional. The surface of the solid is at y = h{x, t) and 
the film is in the y > region with the film-substrate 
interface at y = 0. The system is invariant in the z 
direction and all quantities are calculated for a section of 
unit width in the z direction. 

We assume that surface diffusion is the dominant mass 



transport mechanism, leading to the following evolution 
of the surface profile |1G] : 
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where F is the free energy of the system and X is a 
positive constant. The free energy is composed of elastic 
and surface terms: 



F = Fei+ dx 7^1 + {dh/dxf 



(2) 



where 7 is the surface tension and F^i is the elastic free 
energy, which includes also any elastic contributions to 
the surface tension. We express F^i as F^i = F,; -|- 
5Fei , where i^^; is the elastic free energy of the zero 
strain reference state, and SF^^i is calculated from linear 
elasticity theory. For each value of x, the reference state 
corresponds locally to a flat film of thickness h{x)\ i.e.. 



F 



(0) 



JdxJ_^'dyfi '{h{x),y), where fi '{h{x),y) is 
the elastic free energy per unit volume of a flat film of 
thickness h{x). 

Due to the lateral variations in this reference state, the 
reference stress does not satisfy the condition of mechan- 
ical equilibrium. However, the necessary corrections van- 
ish in the limit a/A — > 0, where a is the length scale over 
which stress varies in the y direction and A is the lateral 
length of typical surface structures. This is because in 
this limit there are no lateral variations in the reference 
stress. As typical experimental islands have A ~ 100 nm, 
and as a is of the order of the lattice constant (see be- 
low) , the corrections to the reference stress are small and 
have been ignored. 

For convenience we work in terms of the reference 
elastic free energy per unit length in the x direction, 



f!,i\H^)) 



S-^' dyfv{h{x),y), instead of the free en- 



ergy per unit volume. As discussed below, the depen- 
dence of f^i on /i is a nonlinear phenomenon and cannot 
be calculated from linear elasticity theory. 

For small strains, the stress is linear in the strain, i.e., 
^ij ~ ^Yj ~^ CijkiCki, where repeated indices are summed 
over. Here Cj,- is the strain tensor, aij is the total stress 



tensor, s^ is the stress in the zero strain reference state 
due to the lattice mismatch and cijki are the elastic coef- 
ficients of the material. According to linear elasticity the- 



ory, dF,i ^ J dx j!^l^ 

of the stress tensor, we find 
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where we have used the inverted Hooke's law Cij — 
SijkiSki- Sijki are the compliance coefficients of the ma- 
terial. Combining Eqs. (0) and d3), we arrive at an ex- 
pression for 5F/6h at the surface: 
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where n is the atomic area of the solid, k is the surface 
curvature, j{9) — j{d) + d^j/dO^ is the surface stiffness 
and 6 is the angle between the normal to the surface 
and the y direction. As Eq. (^ gives SF/Sh at the solid 
surface, all variables in the equation are also given at 
the surface. Both df^^ /dh and s^" must vanish when 
h < 0, since then the film is absent. In principle, Eq. 
(Q) should also contain derivatives of 7 with respect to h. 
However, we believe that the variation of surface tension 
with h away from a step dependence is due to elastic 
effects. Since we included all elastic contributions in the 
zero-strain elastic free energy, we modeled 7 as a step 
function, taking the value of the substrate surface tension 
for h < and the film surface tension for h > 0. Thus 
all partial derivatives of 7 with respect to surface height 
vanish and were omitted from Eq. (0). 

Equations (|^) and (Q) form a complete model of film 
evolution. In order to solve this model one has to eval- 
uate 7(0), Sijki, f^i and s™. The first two are material 
properties, while the last two are properties of the refer- 
ence state, from which one can also calculate the stress 
tensor ctjj using linear elasticity theory. Before estimat- 
ing these quantities we present the results of the linear 
stability analysis of an isotropic flat fllm of thickness C. 
The analysis was carried out using a method similar to 
those used in ||ll| for an infinite film. The height of the 
perturbed film takes the form h{x,t) = C + 6{t)smkx. 
We assumed that the force on the surface due to sur- 
face tension is negligible compared to the force due to 
mismatch stress, and that the stress atj vanishes deep 
in the substrate. Using linear elasticity theory and Eq. 
(Q) with these assumptions, we calculated 5F/6h to first 
order in the perturbation, and combined the results with 
the general evolution equation (||) to obtain the following 
equation for the evolution of 5{t): 
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J h=C 



where 70 = j{9 = 0), M is the plain strain modulus de- 
rived from the elastic constants of the isotropic material, 
and r]{h) is s^ at the surface of a flat film of thickness 
h. s^y vanishes because the flat film is hydrostatically 
strained, and s"i — since in the reference state the 
force on the surface in the y direction vanishes. 

Equation M) implies that the flat film is stable at all 
perturbation wavelengths as long as 
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and the equality holds at the critical wetting layer thick- 
ness. 7o is positive if = is a surface seen in the equilib- 
rium free crystal |12]. At a perfect facet, 70 -^ 00. Hence, 
as 6* = is a facet direction for most of the materials used 
in epitaxial films, 70 is large and positive. Therefore, a 
linearly stable wetting layer of finite thickness can exist 
only if d^y^i /dh? > 0. Note that 77 depends linearly 
on the lattice mismatch e, and hence the l.h.s. of (ph is 
proportional to e'^ , while the r.h.s. of (ph is proportional 
to e^ due to the dependence of /^^ ' on lattice mismatch. 
Therefore, if d'^fj^i /dh'^ > 0, the thickness of the wet- 
ting layer increases with decreasing lattice mismatch and 
diverges in the limit e — > 0. 

Having recognized the importance of the elastic free 
energy of the reference state, /^j , and its dependence 
on film thickness, we now turn to estimate it. This free 
energy depends strongly on the mismatch stress s^, and 
its dependence on the y coordinate. As a result of the 
sharp interface between the substrate and the film, we 
expect s^ to behave as a step function of y with small 
corrections due to elastic relaxation. If we ignore these 
small corrections, the resulting free energy /^^ , is pro- 
portional to film thickness, and its second derivative van- 
ishes. Hence, according to Eq. (^, the thickness of the 
critical wetting layer vanishes. The correction due to 
elastic relaxation is therefore extremely important. It 
turns out that this correction vanishes within linear elas- 
ticity theory. This led some investigators Q to claim that 
the variation in free energy over the interface was due 
to non-elastic effects, e.g. film/substrate material mix- 
ing over the interface. However, we claim that this is 
not necessary, since nonlinear elasticity can explain the 
corrections to the step-function form of the free energy. 

Ideally, first principles, substance-specific calculations 
should be performed in order to evaluate ?/(/i) and 
f^i (h), and we intend to carry out such calculations in 
the future. However, the qualitative general behavior of 
/g; (h) can be obtained from much simpler models. To 
demonstrate this point we carried out the calculation for 



two dimensional networks of balls and springs of vary- 
ing lattice type and spring constants. In this model the 
balls are connected by springs which obey Hooke's law. 
The natural spring length had a step variation over the 
interface. The film underwent a hydrostatic transforma- 
tion strain so that its lattice constant became that of the 
substrate. The network was then allowed to relax whilst 
being constrained in the x direction and free in the y di- 
rection, so that the system boundaries in the x direction 
were fixed to the natural substrate length. 

We calculated the mismatch stress within the film and 
at the film surface for films of varying thickness. How- 
ever, we decided to use the step function form of mis- 
match stress, ri{h > 0) = Me, where Me is the mismatch 
stress in an infinite film, as variations in ?/ only slightly 
altered the wetting layer thickness predicted from (g). 
We also calculated the nonlinear elastic free energy of 
the relaxed system per unit length in the x direction for 
various film thicknesses. A typical behavior of d/ei /dh 

is shown in Fig. I, where it is seen that f^i (h) indeed 
depends on the thickness h. Moreover, the model pre- 
dicts that cP f^J^/dh? > 0, and therefore according to 
the inequality (g|) and the discussion following it, there 
should be a linearly stable wetting layer, whose thickness 
is finite and increases with decreasing lattice mismatch. 
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FIG. 1. Variation with film thickness of the elastic free en- 
ergy of a relaxed balls and springs system, df^^ /dh, as a 
function of film thickness h. The free energy is normalized to 
the infinite film linear elastic energy density, iMe^. hmi is 
the thickness of one monolayer. 

The dependence of /^^ on h arises within our model 
from the elastic relaxation at the surface and its coupling 
to the relaxation at the interface between the substrate 
and the film. A similar effect should occur in real systems 
due to surface reconstruction, for example. 

Our calculations indicate that the general qualitative 
behavior of the mismatch stress and the elastic energy are 
not sensitive to the lattice structures and the values of the 
spring constants. While the detailed behavior close to the 



substrate-film interface (^ 3 monolayers) varied between 
different networks, it showed the same general behavior. 
In all systems (Pf),i/dh? showed exponential decay with 
a decay length of about a monolayer away from the inter- 
face. For the calculations used later in this paper we used 
the function df^°^/dh = Me^ [l _ 0.05exp(-V/im/)] /2 
for /i > and df^^ /dh = for h < 0. hmi is the thick- 
ness of one monolayer. The factor Afe^/2 on the r.h.s. 
ensures that the above expression becomes exact for a 
film of infinite thickness. 

Combining this behavior of df^i /dh with the inequal- 
ity (ph, we obtained an expression for the linear stability 
wetting layer thickness, h^'. 



he /hmi = max {l, ;n [yQ/{AOMe'^hmi)] } 



(7) 



Thus, the wetting layer thickness increases with decreas- 
ing lattice mismatch, as observed in experiments. 

In previous works |^,||J|] on the physics of the wetting 
layer it was assumed that the reference state energy vari- 
ation is a smooth function of h, mainly in order to avoid 
non-analyticities at the interface. In contrast, our ref- 
erence state energy variation behaves as a step function 
of the surface height with a small correction. We have 
shown that the non-analytic behavior at the interface is 
realistic and that the smooth elastic energies in [§|,||j|] 
are unphysical. Tersoff [|| in effect calculated df^i /dh 
via an effective atomic potential for Si/Ge under 4% lat- 
tice mismatch. However, he did not address the stability 
of the flat film to small perturbations. Nevertheless, as 
the function he obtained had d^f^^ /dh^ > 0, a positive 
wetting layer thickness can be predicted from his results. 

In order to model the early evolution of faceted is- 
lands, and to study the effect of an anisotropic form 
of surface tension on the wetting layer, we used the 
cusped form of surface tension given by Bonzel and 
Preuss iQ, which shows faceting in a free crystal: ^{6) = 
7o[l-H/3|sin(7r6i/(26lo))|], where /3 « 0.05 and 6*0 is the 
angle of maximum 7. The value of 70 was taken as 1 
J/m^ in the substrate and about 75% of that in the film 
(as is the case for Si/Ge). This ensures a wetting layer 
of at least one monolayer. We considered a crystal which 
facets at 0°,±45° and ±90° with 6*0 = tt/S. The cusp 
gives rise to 7 = 00 and hence all faceted surfaces will 
have an infinite linearly stable wetting layer. However, a 
slight miscut of the low-index surface leads to a rounding 
of the cusp, which can be described by 



7^=70 1 + /3a sin^( 



'200 



G- 



(8) 



where, for example, G — 500 corresponds to a miscut 
angle, A9 « 0.1°. 

According to Eq. (R), anisotropic surface tension 
greatly enlarges the linearly stable wetting layer thick- 
ness. Does this conclusion survive beyond linear stabil- 
ity analysis? When a linearly stable flat film is perturbed 



strongly so that the surface orientation in some regions 
is far from the 9 = direction, the local surface stiffness 
in these regions is much smaller than the 6 = stiffness. 
This tends to destabilize the linearly stable film. Indeed, 
we carried out Monte Carlo simulations that showed that 
films thinner than the linear wetting layer were unstable 
to random perturbations greater than a certain critical 
amplitude (see Fig. 2). The linear elastic energy was 
calculated by the method used by Spencer and Meiron 
WM. Hence films thinner than the linear wetting layer 
thickness are metastable. When large perturbations were 
applied, faceted islands developed in the film, which un- 
derwent Ostwald ripening at later stages of the evolution. 
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FIG. 2. Evolution of a randomly perturbed film, in which 
perturbations were larger than the critical perturbation am- 
plitude. Lattice mismatch in this film is 4%. The initial film 
surface is shown as a thin solid line. The dashed line shows the 
film surface at a later time. The linear wetting layer thickness 
is shown as a thick solid line. 

The critical perturbation amplitude was found to be 
proportional to £~^. It was largely independent of film 
thickness for h < he as long as perturbations did not 
penetrate the substrate, when it became much larger. It 
was also largely independent of cusp smoothness G, un- 
like the linear wetting layer thickness which depended 
strongly on G. The size of the critical perturbation am- 
plitude in monolayers is plotted as a function of lattice 
mismatch in Fig. 3. The linear wetting layer thickness 
for G = 500, M = 1.5 x lO^A^/m^ and h„u = 5A is also 
shown for comparison. When the lattice mismatch is 
small, the critical perturbation amplitude is much larger 
than a monolayer. Hence, in practice, flat films thinner 
than the linear critical thickness are stable at small lattice 
mismatch. On the other hand, for large mismatch a per- 
turbation smaller than a monolayer is sufficient in order 
to destabilize the linearly stable wetting layer. There- 
fore, in practice, the wetting layer will be a single mono- 
layer at large lattice mismatch. Our predictions cannot 
be compared with current experiments involving the wet- 
ting layer and lattice mismatch variation [E|, since they 



were carried out with deposition flux. We hope this work 
will encourage such experiments to be performed, and we 
are currently adding deposition to our model. 
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FIG. 3. Variation of critical perturbation amplitude with 
lattice mismatch 
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